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Ferromagnetic couplings in spin clusters are shown to 
be strongly enhanced compared to those for an ordered im- 
purity arrangement, even for the same spin separation and 
hole doping. The consequent energy-enhancement of the 
cluster-localized spin-wave modes indicates a potentially sig- 
nificant role of positional Mn disorder in enhancing Tq. 
Within a simple model involving two spin-excitation energy 
scales corresponding to weakly and strongly coupled spins, 
the temperature dependence of magnetization is found to be 
in good agreement with the SQUID magnetization data for 
Gai_a;Mna;As. 



I. INTRODUCTION 

Following the discovery of ferromagnetism in Mn- 
doped III-V semiconductors such as p-type Ini_j;Mnj;As 
and Gai_a;Mna;As, there has been considerable interest 
recently in understanding the nature of the ferromagnetic 
state and the fundamental mechanism of spin coupling 
in these dilute magnetic semiconductors (DMS). With a 
highest transition temperature (Tc) of llOK for Mn con- 
centration X = 0.053, [1] and a realistic possibility of 
room temperature ferromagnetism in Gai-^^Mn^^N, with 
a highest reported Tc value of 940K, [2] these materi- 
als have potential applications in seamlessly integrated 
non- volatile semiconductor memories. 

Mn doping in III-V semiconductors plays a dual role 
of providing magnetic ions as well as acceptor sites, and 
the ferromagnetic interaction between the 5 = 5/2 Mn++ 
ions is mediated by the itinerant charge carriers, which 
are antiferromagnetically coupled to the Mn spins. The 
hole doping is heavily compensated by As antisite defects, 
resulting in a carrier concentration much smaller than 
the magnetic impurity concenteration [p << x), no that 
the DMS system provides an interesting compliment to 
the Kondo system. In addition, the involvement of low- 
momentum, long-wavelength carrier states at the top of 
the valence band in the spin coupling process results in 
physical simplification and qualitative independence on 
details of the electronic band structure. 

Various approaches have been employed to study the 
carrier-mediated ferromagnetism in DMS. These include 
the mean- field (Zener model), [3-8] RKKY interaction, 
[1] spin-wave theory, [9,10] dynamical mean field the- 
ory, [11] Monte Carlo simulations, [12,13] a generalized 
RKKY approach which takes into account the spatial 
variation of the impurity-induced carrier spin polariza- 
tion beyond linear response, [14] and a mean-field-plus- 
fiuctuation (MF-I-SF) approach within a Hubbard-?/ rep- 
resentation of the magnetic impurities. [14] 



Recently, there has been interest in understanding the 
role of the positional Mn disorder on macroscopic prop- 
erties such as temperature dependence of magnetization, 
susceptibility, specific heat etc. [15-19] Relevance of clus- 
tered states in dynamical and transport properties has 
also been recently studied. [20] In this paper we study the 
distribution of ferromagnetic spin couplings arising from 
clustering of Mn spins, and the consequences on collec- 
tive spin-wave excitations which play an important role 
in the low-temperature thermodynamics of magnetic sys- 
tems. We use the numerical MF-I-SF approach in which 
the Mn disorder is treated exactly and electron correla- 
tion effects are teated within the random phase approxi- 
mation. 



II. HUBBARD-t/ REPRESENTATION 

We consider a purely fermionic (Hubbard- [/) represen- 
tation for the randomly distributed magnetic impurities 
on a cubic host lattice: 



<ij>cr 



<Ij>a 



Here i,j refer to the host sites, / to the magnetic impurity 
sites, Cd is the impurity on-site energy and nj — {hj^ + 
fill) / 2 is ti^s spin-averaged impurity charge density. For 
simplicity, we take the same hopping [t' = t) between 
the host-host and host-impurity nearest-neighbour pairs 
of sites. The energy-scale origin is set so that the host 
on-site energy is zero. We take the impurity level to lie at 
the top of the host band (e^ = 6), and U — 4 throughout. 

In the Hartree-Fock (mean-field) approximation, the 
interaction term reduces to a magnetic coupling of the 
electron to the local mean field A/: 



HF 



(2) 



4'|[ct]^'/ in terms 



— (a||a||), and the mean field A/ is 



where the electronic spin operator aj 
of the spinor — 

self-consistently determined from the ground-state ex- 
pectation value 2A/ = [/((?/). Thus, in the classical 
(Hartree-Fock) limit, the interaction term reduces to the 
corresponding form of the double-exchange interaction 
— JS'/.f?/, with the mean field A/ — J{Si) representing 
the impurity-induced local magnetic field. 
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FIG. 1. Variation of impurity magnetization on selected 
impurity sites with iterations, showing crossover to an inho- 
mogeneous ferromagnetic state with vanishing magnetic mo- 
ment on a fraction of impurity sites. 

We consider a cubic host lattice (spacing a = 1) 
with N = 8^ sites, and focus on a disordered impu- 
rity arrangement with Ni^p = 27 magnetic impurities 
{x w 5%) placed at locations 1,4,6 in all three direc- 
tions. This arrangement results in an isolated impurity 
spin at location (1,1,1) with nearest-neighbour separa- 
tion of 3, and an impurity cluster of 8 spins at locations 
(4,4,4),(4,6,4),(6,4,4),(6,6,4),(4,4,6),(4,6,6),(6,4,6),(6,6,6) 
with nearest-neighbour separation of 2. For comparison, 
we also consider an ordered impurity arrangement of 64 
impurities (x = 1/8) placed on a cubic superlattice (spac- 
ing ttsup = 2). 

The undoped (insulating) state corresponds to elec- 
tron fillings = N and = N — Ni^^; all minority- 
spin (I) impurity states (pushed out of the host band by 
Coulomb repulsion) are then unoccupied at T = 0, re- 
sulting in local-moment formation on all impurity sites. 
Hole doping is introduced by reducing JV| , and band fill- 
ings are so chosen that the Fermi energy lies in gaps 
between (nearly) degenerate groups of eigenvalues. We 
take Ni = 502, corresponding to p « 2%. 

The variation of impurity magnetization in the self- 
consistency process (Fig. 1) reveals an interesting slow 
dynamics associated with longitudinal fluctuations. A 
nearly homogeneous ferromagnetic state is obtained ini- 
tially which abruptly becomes unstable (around the 10**^ 
iteration) towards an inhomogeneous state with vanish- 
ing impurity moment on a finite fraction of impurity 
sites. Interestingly, these effectively non-magnetic impu- 
rity sites are neighbours of the isolated spin at location 
(1,1,1). At the 20**^ iteration, this inhomogeneous ferro- 
magnetic state is nearly self-consistent, and we take it as 
the HF state for the calculations. The stability analysis 
[14] yields fAmax = 1-0007, confirming the near stability. 
The above behaviour indicates an unusual susceptibility 
of the DMS towards large-amplitude, slow fluctuations in 
the impurity magnetic moment. 

We now consider the consequences of the spin cluster- 
ing, starting with the impurity states. For a single mag- 
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FIG. 2. Magnetization on impurity sites, showing a 
distinct reduction on the strongly coupled cluster sites 
(14,15,17,18,23,24,26,27). 

netic impurity, the minority-spin (J,) impurity hole state 
(pushed out of the band) is exponentially localized. 
For multiple impurities, the overlap between impurity- 
state wavefunctions of neighbouring impurities leads to 
a mixing which can be represented by an effective hop- 
ping t' between impurity sites. For the ordered impu- 
rity arrangement, all neighbouring hopping terms and 
impurity-state energies are identical, and resonant hop- 
ping results in impurity-band formation. For the disor- 
dered case, the large variation in the effective impurity 
hopping t' between cluster sites and isolated sites leads 
to localization of impurity states within clusters. 

As the mixing between neighbouring impurity states 
is mediated through the intervening host sites, the hole 
density hi is somewhat enhanced on the cluster host 
sites. The remaining hi on the magnetic sites is there- 
fore somewhat reduced, and the corresponding increase 
in the electron density rtj = 1 — hi implies a lowered 
magnetization of cluster spins. This picture is confirmed 
in Fig. 2, showing that the (HF) impurity magnetization 
mj ~ nil ~ cluster spins is distinctly lower than 

that of isolated spins. The lowered hole density hi on 
cluster sites also reduces the diagonal terms [x°]// of the 
zeroth-order, particle-hole propagator [Eq. (4)], indicat- 
ing stronger spin couplings. 



III. FERROMAGNETIC SPIN COUPLINGS 

The ferromagnetic spin coupling between impurity 
spins at sites / and J can be obtained as [21] 



i/2[x°(a; = 0)] 



ij , 



(3) 



where [x'^{uj)]jj is the zeroth-order, particle- hole prop- 
agator, evaluated in the self-consistent ferromagnetic 
state: 
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FIG. 3. The spin couplings Jij = U'^[x°]iJ for all nearest 
neighbours of impurity spins, showing a strong enhancement 
in the cluster couplings compared to the spin coupling for a 
ordered impurity arrangement (dashed line). 

The spin couplings are shown in Fig. 3 for all six near- 
est neighbours of each impurity spin, the separations be- 
ing either 2 or 3. Interestingly, the spin couplings are 
grouped into three distinct classes — weak, intermedi- 
ate, and strong. A small fraction of the couplings arc also 
weakly negative (antiferromagnetic), indicating compet- 
ing interactions. Also shown is the NN spin coupling 
for the ordered impurity arrangement (A^imp = 64), ex- 
trapolated to the same hole concentration. The dramatic 
enhancement of cluster spin couplings, although the NN 
separation is the same (2), is beyond the simple RKKY 
picture wherein the spin coupling J^Xu depends only on 
the carrier concentration and spin separation. 



IV. SPIN- WAVE EXCITATIONS 

The (time-ordered) spin-wave propagator 

X7/{t - t') = {^o\T[Si{t)S+{t'mo) , (5) 

involving the spin- lowering {S^) and spin-raising {Sj) 
operators at magnetic impurity sites I and J, describes 
the low-energy transverse spin fluctuations about the 
HF state. In the random phase approximation (RPA), 
[x '"(w)] = [x^{u))]/l — U[x^{u))], and we obtain the spin- 
wave energies from the pole condition 1 — f/A„(tJ„) = 0, 
where A„(w) are the eigenvalues of the [x*'(a;)] matrix, 
evaluated in the self-consistent state. 

The localization-induced enhancement of NN spin cou- 
plings in the disordered impurity configuration also leads 
to substantially higher spin-wave energies. Figure 4 
shows a comparison of the spin- wave energies a;„ for the 
disordered and ordered impurity arrangements at the 
same hole concentration (p « 2%). For the disordered 
arrangement, three localized modes associated with the 
three non-magnetic impurity sites have been excluded. 
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FIG. 4. Strong stiffening of the high-energy spin-wave 
modes due to clustering. 
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FIG. 5. Variation of the Participation Ratio for the 

spin-wave modes with their energy, showing that the stiffening 
of the high-energy spin- wave modes is due to their localization 
over the strongly-coupled cluster spins. 

For the ordered arrangement, energies for the nearest 
filling (A^i = 493, p = 3.7%) have been (linearly) ex- 
trapolated to p = 2%. The nearly uniform energy in- 
terval in this case indicates dominant nearest-neighbour 
spin coupling, for which an energy interval of zJS/3 fol- 
lows from the spin- wave energies Wq = zJS{l — 7q), 
corresponding to plane-wave modes of momentum q = 
(2tt/ L){nx,ny,nz), with L = A. It is clear that while the 
low-energy modes are softened, the high-energy modes 
are significantly stiffened by impurity disorder. 

To examine the spatial nature of the spin-wave modes 
(wave function (/)^), we study their Participation Ratio 
Pn = [E/('/'n)^]V[E/(</'n)'']> which provides a quanti- 
tative measure of the number of sites over which the 
wavefunction is spread. Figure 5 shows that the stiffened 
high-energy modes are indeed cluster-localized {PR ^ 8). 
On the other hand, the low-energy modes show a large 
range of PR values, indicating states of both extended 
and localized nature, suggesting different mechanisms of 
softening, as discussed below. 
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FIG. 6. Temperaturo-dcpondonce of M{T) from Eq. (7), 
along with SQUID data for Gai-a:MnxAs (circles). 

We now consider the temperature dependence of the 
magnetization in a DMS having a finite fraction (/i) of 
weakly coupled spins, whose thermal disordering satu- 
rates at a relatively lower tcmperatiire, resulting in a 
concave M(T) behaviour. While a full renormalized spin- 
wave-theory evaluation of M{T) is possible, we consider 
a simpler model in which the spin excitations arc char- 
acterised by two energy scales uj\ and 0)2, corresponding 
to the weak and strong couplings, respectively. Taking 
mean- field values as 1, the reduction in magnetizations 
Ml and M2 of weakly and strongly coupled spins, due to 
the bosonic spin excitations, are obtained as: 



Ml = 1 - 
M2 = 1 



Ml 



exp(M2a;i/fcBT) - 1 
M2 



CXp(M2W2/fcBT) - 1 ' 

yielding the average magnetization 

M(r) = /iMi + (1 - /i)M2 . 



(6) 



(7) 



Here the renormalized excitation energies M2iiJi and 
M2a;2 reflect the role of the strongly coupled spins in the 
overall magnetic ordering. The average magnetization 
M(T) is shown in Fig. 6, along with the SQUID magne- 
tization data for Gai-a^Mn^jAs with x = 0.053, [22] with 
best-fit parameters /i = 0.55, oji = 15K, and W2 = 120K. 



V. CONCLUSIONS 

The positional Mn disorder in DMS results in iso- 
lated spins and spin clusters, with substantially differcnit 
impurity-state overlap and effective impurity-impurity 
hopping. This large hopping disorder leads to localiza- 
tion of impurity states over different clusters. The fer- 
romagnetic spin couplings J/j = U^X^j of isolated 
impurity spin are consequently appreciably weakened, 
whereas the enhanced overlap of the cluster-localized im- 
purity states significantly increases the cluster spin cou- 
plings. The consequent disorder in the [x°] matrix leads 



to localization of spin waves over different clusters. The 
low-energy spin-wave modes arc softened due to i) lo- 
calization on the weakly-coupled isolated spins, or ii) ef- 
fect of disorder-induced competing AF interaction on the 
long-wavelength modes. On the other hand, the local- 
ization of high-energy modes over strongly coupled spin 
clusters results in a significant enhancement of their en- 
ergy. The temperature dependence of magnetization, us- 
ing two spin-excitation energy scales corresponding to the 
weakly and strongly coupled spins, is in good agreement 
with the SQUID magnetization data for Gai-xMn^As. 
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